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Abstract 

^ | We discuss a non-perturbative lattice calculation of the ghost and gluon propaga- 

tors in the pure Yang-Mills theory in Landau gauge. The ultraviolet behaviour is 
checked up to NNNLO yielding the value = 269(5)+^ 2 MeV, and we show 



o 



that lattice Green functions satisfy the complete Schwinger-Dyson equation for the 
ghost propagator for all considered momenta. The study of the above propagators 
at small momenta showed that the infrared divergence of the ghost propagator is 
enhanced, whereas the gluon propagator seem to remain finite and non-zero. The 
result for the ghost propagator is consistent with the analysis of the Slavnov- Taylor 
O ■ identity, whereas, according to this analysis, the gluon propagator should diverge 

Q , in the infrared, a result at odds with other approaches. 

CPHT-PC076.1205, LPT-06-02, UHU-FT-18 

X 

1 Introduction 

In this talk we report on a non-perturbative lattice study of basic correlation functions of 
the Euclidean Landau gauge pure Yang-Mills gauge theories, i.e. the gluon propagator 

GT^-P) = ^{^- P -f)^ (!) 
and the ghost propagator 

F (2)ab^_ p)=S a b FiP)_. ( 2 ) 
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These propagators have been successfully checked by perturbation theory up to NNNLO 
in the ultraviolet domain, (PQ,[2]) yielding the non-perturbative value of 

A^ =0 = 269(5)^ 2 MeV. (3) 

Now we concentrate on the infrared exponents ap and ag that describe power- law devi- 
ations from free propagators when p —>■ 

G(p 2 ) oc (p 2 ) aG F(p 2 ) oc (p 2 ) ap . (4) 



2 Non-perturbative study of the infrared exponents 

OL F , ®G 

Diverse analytical approaches (study of truncated Schwinger-Dyson equations and of 
renormalisation group equation [Hj, [I]) agree that the infrared divergence of the ghost 
propagator is enhanced, i.e. op < 0; while they predict different values for a G , mostly 
around ~ 1.2. This means that the gluon propagator is suppressed in the infrared. 
Lattice simulations confirmed the prediction for the ghost propagator, whereas the lattice 
gluon propagator seems to remain finite and non-zero in the infrared, i.e. «g = 1 [HJ • 

Let us consider the Slavnov- Taylor identity [Hj relating the three-gluon vertex Tx^, 
the ghost-gluon vertex T\^(p, q; r) 3 and the propagators: 

p X T x ^(p,q,r) = ^3y(^r 2 - r A r„)f A/i (r, p; q) - ^j^y(<^ 2 - 1x1$ \»{<l,V\ r). (5) 

Taking the limit r — > keeping q and p finite, and using G{r 2 ) ~ (r 2 ) a ° , one finds the 
following limits on the infrared exponents [2] 

«g < 1 gluon propagator diverges in the infrared, and , . 

«f < the divergence of the ghost propagator is enhanced in the infrared. 

The limit on disagrees with many other analytical predictions [3], except some cases 
in Let us try to understand this discrepancy. All these methods rely on a commonly 
accepted relation between the infrared exponents 

2a F + a G = 0. (7) 

which we shall discus now. The origin of this relation is the dimensional analysis of the 
Schwinger-Dyson equation for the ghost propagator (SD): 



1 2at f d 4 q ( F(q 2 )G{{q-kf 



l + g$N c 



(k ■ g) 2 - k 2 q 2 
k 2 (q — k) 2 



Hi(q,k) 



F(k) " J (2tt) 4 y q 2 (q-k) 2 

where Hi(q, k) is one of the scalar functions defining the ghost-gluon vertex: 

q v '^v'u{-q, k;q- k) = q^H^q, k) + (q - k) v H 2 (q, k). (9) 
3 r is the momentum of the gluon, q is the momentum of the entering ghost. 
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The large momentum behaviour ([S],jS]) of this vertex depends on the kinematic configu- 
ration: 

^•rjf(-p,0;p) = l to all orders 

^•r^ s (-p,p;0) = i + ^(P 2 ) + ... 1 j 

Note that in the case of the vanishing momentum of the out-going ghost (and only in this 
case) one has 

H 1 (q,0)+H 2 (q,0) = l. (11) 

If both H 12 are non-singular then one can suppose Hi(q, k) ~ 1 in (JSJ), and (J7J) is straight- 
forward by a dimensional analysis. However, we have a priori no reason to think that the 
scalar functions Hi(q, k) and H2(q, k) are separately non-singular for all q, k. Writing, for 

example, Hi(q, k) ~ {(f) ar hi \^r, ^j, with a non-singular function hi, we keep all the 
generality of the argument admitting a singular behaviour of the scalar factor Hi(q,k). 
Doing the dimensional analysis of eq. (JEJ) without putting Hi(q,k) ~ 1, we obtain that 
the relation (JZJ) is satisfied if and only if the following triple condition is verified JZJ: 



2a F + a G = { a r = (12) 




The first and the last conditions are compatible with limits coming from the analysis of 
the Slavnon- Taylor identity ©, and are also consistent with our lattice simulations [7j. 
If one of the conditions (fT2|) is not verified then (JZj) should be replaced by 

2ap + ac + «r = 0. (13) 

In the following section we present the results of a numerical test of the relation and 
thus we probe the validity of the condition on ar- 



3 Lattice check of analytical predictions 

As we have already mentioned, lattice Green functions are checked at high precision with 
perturbation theory at large momentum. In order to test the validity of lattice predictions 
at small momenta we verified that lattice Green functions satisfy the ghost SD equation 
written in the form 

F(P 2 ) = 1 + 90^/^(^(0) ■ Fffi(Ap)>, (14) 

where F^ ^(A, p) is the correlator of the ghost and anti-ghost fields in a background field 
A. Note that the three-point function in the r.h.s is in a mixed coordinate-momentum 
representation, and that here we make no assumption about the vertex. The traditional 
form (jHJ differs from ()14j) by a Legendre transformation. We see from Fig^ that the 
lattice Green functions match pretty well the SD equation (j!4)l in both the ultraviolet 
and infrared regions. Now we try to check the validity of the non-singularity assumption 
for the scalar factor Hi(q, k). For this purpose we check numerically whether our lattice 
propagators satisfy the equation (JSJ) completed with the assumption Hi(q,k) = 1 (cf. 
Fig|21 [7J). We see that at small momenta (below m 3 GeV) the SD with the assumption 
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Figure 1: Checking that lattice Green functions satisfy the ghost SD equation (|14|l . The l.h.s vs r.h.s 
of (Q3|) is plotted at the left, and at the right wc plot F(p 2 ) - goT^f a ^(A°(0).F£££(A,p)) compared 
to 1 




Figure 2: Checking whether lattice Green functions satisfy the ghost SD equation JHJ with an assumption 
Hi(q, k) = 1. The upper line(circles) correspond to the loop integral in (JSJ, and the down line (triangles) 
corresponds to 1/F(p 2 ) — 1. In this plot a -1 w 3.6 GeV. 



H\(q,k) = 1 is not satisfied. This suggests that the scalar function H\(q,k) plays an 
important role in the infrared gluodynamics. Finally, a direct test of the relation (JJJ) and 
thus, as we have seen in f)12jl . of the value of ar can be done. For this we plot the quantity 
F 2 {p 2 )G{p 2 ) FigO If all the conditions (fT2^) are satisfied this quantity should be constant 
in the infrared. We see from FigJSlthat in the infrared (below w 600 MeV) the quantity 
F 2 G is not constant, and thus one of the conditions (jl2j) is not verified. We have seen 
that the conditions c\p 7^ and ap + ac < 1 are consistent with the limits © from the 
Slavnov- Taylor identity (JHJ). We have also seen (cf. Fig ll|2jl that neglecting the vertex 
is not possible in the infrared, because in this case the ghost SD equation is no longer 
satisfied by lattice propagators. Thus the only possibility is to admit that H 1 (q, k) plays 
an important role, and that the relation (J7J) is not verified. The modified form (|13|) that 
takes in account the singularity of H\{q,k) should be considered (according to FigEJ), 
with ar < in our parametrisation. 
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Figure 3: Direct test of the relation 2cxf + olq = 0. If the last is true F 2 G has to be constant in the 
infrared. We see that it is clearly not the case. In these plots a -1 « 1.2 GeV, so the peak is located at 
w 600 MeV. 

4 Conclusions 

We have seen that lattice simulations allow a whole momentum range study of Green 
functions of a non-Abelian gauge theory in Landau gauge. They have been tested by per- 
turbation theory up to NNNLO at large momentum. We have also checked (numerically) 
that lattice Green functions satisfy the complete ghost Schwinger-Dyson equation (JT4"jl for 
all considered momenta. These tests allow us to conclude that numerical simulation on 
the lattice give relevant results not only in the ultraviolet domain but also in the infrared 
one. 

Our analysis of the Slavnov- Taylor identity showed that the power-law infrared diver- 
gence of the ghost propagator is enhanced in the infrared (compared to the free case), and 
that the gluon propagator must diverge in the infrared. The latter limit is in conflict with 
most present analytical estimations j^j, that support a vanishing gluon propagator in the 
infrared. In [4] the author replaced H\ (q, k) by its perturbative expansion, and an interval 
of values of «g was found, including those compatible with (jHJ. Lattice simulations point 
to a finite non-vanishing infrared gluon propagator, in conflict with the above analysis of 
the Slavnov- Taylor identity. 

Our numerical studies showed that the commonly accepted relation (J7J) between the 
infrared exponents is not valid, because F 2 G is infrared suppressed, and hence 2otF + ctG > 
0. This statement is supported by the fact that lattice propagators do not match the 
reduced SD equation (FigEJ), whereas the complete one is perfectly verified (FigHJ). This 
speaks in favour of a singularity in the scalar factor Hi. Note that it does not contradict 
the non-renormalisation theorem jB] which implies for the renormalisation constant of the 
ghost-gluon vertex : = 1. 
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